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A transmission problem across a fractal self-similar interface 

Yves Achdou * Thibaut Deheuvels^ 


Abstract 


We consider a transmission problem in which the interior domain has infinitely ramified structures. 
Transmission between the interior and exterior domains occurs only at the fractal component of the 
interface between the interior and exterior domains. We also consider the sequence of the transmission 
problems in which the interior domain is obtained by stopping the self-similar construction after a 
finite number of steps; the transmission condition is then posed on a prefractal approximation of the 
fractal interface. We prove the convergence in the sense of Mosco of the energy forms associated with 
these problems to the energy form of the limit problem. In particular, this implies the convergence 
of the solutions of the approximated problems to the solution of the problem with fractal interface. 
The proof relies in particular on an extension property. 

Emphasis is put on the geometry of the ramified domain. The convergence result is obtained when 
the fractal interface has no self-contact, and in a particular geometry with self-contacts, for which an 
extension result is proved. 


1 Introduction 


Transmission problems naturally arise in various fields of physics and have been extensively studied. An 
introduction to this class of problems can be found in [20] and several applications are detailed in [14] . 
Such problems were more recently studied in the case when the interface is irregular, Lipschitz contin¬ 
uous or even fractal. These problems find many applications, such as the study of rough electrodes in 
electrochemistry or diffusion across irregular membranes in physiological processes, etc. (see Ell USD- 
Several transmission problems with fractal interfaces have been studied in the case of the Koch flake or 
the Sierpinski gasket in 2D and 3D (see e.g. [251 EH ES3 ESI [HI EH IIS USD- 

This paper deals with transmission problems between two domains flj n t and D ext where is a 
ramified bounded domain as defined in Section 12721 (see Figure [2]). The domain Dj nt presents infinitely 
many ramifications, and its boundary contains a fractal self-similar set T which plays the role of the 
interface. 

The domain Di nt can be seen as a bidimensional idealization of the bronchial tree, for example. Since 
the exchanges between the lungs and the circulatory system take place only in the last generations of the 
bronchial tree (the smallest structures), it is reasonable to consider transmission problems with specific 
transmission conditions accross the fractal boundary T of the tree. We will however limit ourselves to 
simple transmission conditions. 

The fractal boundary T belongs to a family of self-similar sets introduced by Mandelbrot et al. in [30l . 
Elliptic boundary value problems in the domain have been studied in [3], and traces and extension 
results for these domains have been proved in mm- 
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The considered problem can be formally stated as 

'-A u = f 

M = o 

[d n u\ = au 
9 n U int = <5n^ext = 0 
dn.Uext — 0, Mint — Mo 
, ^n^ex t = 0 


where a is a positive real number, D is a regular bounded open domain in the plane with D = Q; nt U f2 ext 
such that rgD, and T° is a line segment included in the boundary of flint- The sets f!j nt and f2 ext are 
disjoint subdomains of D , and [u] ( resp. [ d n u ]) denotes the jump of u ( resp. of the “normal derivative” 
d n u of u) accross the fractal set T. Since the interface T is fractal, the normal derivative on T has to be 
understood in a suitable weak sense, which will be made precise later. 

Problem (0 is a model problem. Its study is the first step in the modelling of physical transmission 
problems in ramified structures. 

The goal is to study approximations (P„) of problem 13), obtained by stopping the construction of the 
ramified domain Qj nt at step n. The interfaces in the problems ( P n ) are called prefractal approximations 
of the fractal set T. They consist of disjoint finite unions of line segments, which makes problems (P„) 
much simpler than problem (P). A natural question is to understand the asymptotic behavior of the 
problems (P n ) as n —> oo, and in particular to investigate the convergence of the solutions u n of the 
problems (P„) to the solution u of problem (P). 

Remark 1. In contrast with the references f28l [37[ [32 29, \TT\, UM 1731 UHL the boundary value problem 
does not only involve transmission conditions at the interface between flint and f2 ex t- there are also 
homogeneous Neumann conditions on the polygonal part of dfiint/ as a consequence, the traces of u[ n t and 
ite X t on this set do not match a priori. Coping with these discontinuities will be a difficulty in studying 
the convergence of the solutions u n of the problems (P n ) to the solution u of problem (P). 

Remark 2. We have chosen that the source term in (P) appear both in the Dirichlet boundary condition 
on To for Mint and in the Poisson equations in Ui n t and S2 ext ; this is of course completely arbitrary. 

A crucial step in the study of the asymptotic behavior of ( P n ) is the question of extending functions 
defined in the domain flint- More precisely, it is of particular importance that fl; n t should satisfy a 
W 1,p -extension property for some p £ [l,oo], i.e. there should exist a bounded linear operator 

£ : WA’^flint) tT 1 ,p (IR 2 ) 

such that £(u) |n . t = u for all u £ lT 1 ’ p (fli nt ). The domains satisfying this property for all p £ [1, oo] will 
often be referred to as Sobolev extension domains. 

It is well known that every Lipscliitz domain in R n , that is every domain whose boundary is locally 
the graph of a Lipschitz function, is a Sobolev extension domain. Calderon proved the kC 1,p -extension 
property for p £ (1, oo) (see [9]), and Stein extended the result to the cases p = 1 and p = oo (see [42]h 
In [22], Jones generalized this result to the class of (e, <5)-domains, also referred to as Jones domain, 
or locally uniform domains (see [3lj). In dimension two, the definition of (e, 3)-domains is equivalent 
to that of quasi-disks, see [32]. This extension result is almost optimal in the plane, in the sense that 
every plane finitely connected Sobolev extension domain is an (e, <5)-domain, see [22] [32]. The case of an 
unbounded domain in R 2 has been studied in [44] . Extension properties for domains which do not have 
the (e, <5)-property have been studied e.g. in [33, where the authors examine in particular the case of 
domains with cusps. 

When the fractal boundary T of the domains flint studied in this paper has no self-contact, fli„t can 
be proved to be an (e, J)-domain, and and is therefore a Sobolev extension domain. However, in the 
case when the boundary self-intersects, flint does not have the lT 1,p -extension property for all p ^ 1. In 
particular, the extension property does not hold for p = 2, which is the relevant case here, since the vari¬ 
ational formulations of (3 naturally involve the spaces iJ 1 (flint)- In the particular geometry considered 


in flint C flext; 
on r, 
on T, 

on <9f2 in t \ T, 
on T°, 


(■ P ) 
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in Section [5l it will be proved that the transmission condition imposed on T yields a better regularity of 
the trace on T, for functions belonging to the function space arising in the variational formulation. It is 
then possible to deduce an extension result in this case (see Theorem O- 

The main question investigated in this paper is the question of the convergence in the sense of Mosco 
of the energy forms associated with problem (P n ) to the energy form of the problem (P). The notion of 
Mosco-convergence, or M-convergence, was introduced in [33], see also [3B]. It is a stronger convergence 
in the space of forms than T-convergence. In particular, it also implies the convergence of minimizers of 
the energy forms to the minimizer of the limit form. The M-convergence of forms is equivalent to the 
convergence of the resolvent operators associated with the relaxed forms in the strong operator topology 
(see [3611. 

The main results of this paper are Theorems [6] and [9] Theorem [G] is about the convergence of the 
energy forms associated with (P„) in the sense of Mosco to the energy form associated with (0, in the 
case when the fractal interface has no self-contact. The proof uses the extension operator from H l (yi mt ) 
to lf 1 (R 2 ) as a main ingredient. The existence of a continuous extension operator from P 1 (Hi nt ) to 
P 1 (R 2 ), for a particular geometry where the fractal interface self-intersects, is stated in Theorem[9l As 
a consequence, the proof of Theorem [6] can be reproduced in this case, to show the M-convergence of the 
energy forms. 

The paper is organized as follows: the geometry of the interior and exterior domains is detailed in 
Section [2j as well as the prefractal geometry. Section [3] is devoted to the study of the function spaces 
involved in the paper, and emphasis is put on trace and extension results for the domains under study. The 
considered transmission problem is described in Section^ Section POl is devoted to the M-convergence of 
the energy forms associated with the problem with prefractal interface to the energy form of the problem 
with fractal interface, in the case when the boundary of Hint has no self-contact. In Section[5l a particular 
geometry in which the fractal part of the boundary of Hint self-intersects is considered; an extension result 
is proved in this particular case and the M-convergence of the energy forms follows. 


2 The geometry 

2.1 The fractal interface 


2.1.1 Definitions 


Consider four real numbers r,/? i,/ 32,0 such that 1/2 ^ r < l/y/2, /?i > 0, /?2 > 0 and 0 < 9 < 7r/2. Let 
fi, i = 1, 2 be the two similitudes in M 2 given by 


h 



/—/3A cos 9 — X 2 sin 

1^/32 ) r sin 9 + X 2 cos 9 J ’ 


h 



//3i\ / x\ cos 9 + x 2 sin 9 

\p2 J Xi sin 9 + x 2 cos 9 


The two similitudes have the same dilation ratio r and opposite angles ±d. One can obtain /2 by 
composing /i with the symmetry with respect to the vertical axis {x\ = 0}. 

Let T denote the self-similar set associated with the similitudes /i and / 2 , i.e. T is the unique compact 
subset of R 2 such that 

r = / 1 (r)u/ 2 (r). 


It was stated in m (see m for a complete proof) that for any 9, 0 ^ 9 < n/2, there exists a unique 
positive number rg € [1/2, \/\/2[ which only depends on the angle 9 such that 


o if 0 < r < rg, then T is totally disconnected, 
o if r = rg, then T is connected. 


In the following, we will always assume that r ^ rg. 
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Notations For every integer n > 0, we note A n = {1,2}". For a £ A ni we note f a the similitude 
f tJl o ... o f & n ■ We agree to extend the notation to the case n = 0: f a = Id if a £ Aq. We also introduce 
the notation A := \J n>0 A n - 

For a £ A, we note T a = f a (T°), and for every integer n 0, F"= U /-( r °)- 

o-eA. 


2.1.2 Hausdorff dimension of T 


If r ^ Tg, then it can be seen that the open set condition (or Moran condition) holds, see [34] or [25] for 
a definition. The open set condition is satisfied e.g. for the domain Q- lnt defined in ©, if Assumption Q] 
below is satisfied (Theorem [2] proves the existence of such a domain). 

The open set condition implies that the Hausdorff dimension of T is 


d := dim# T = 


log 2 
logr 


see |[T[ .25!.. If 0 ^ 8 < 7r/2, then 1/2 ^ r ^ rg < l/\/2 and thus 1 ^ d < 2. 
In the case when r = r^, introduce the set 


s = /i(F) n/ 2 (r). ( 1 ) 

In this case, the fractal set T self-intersects, and union of the images of S by the similitudes f ai o... o f an , 
<7i,... , a n £ {1,2} is the set of the multiple points of T. 

Two situations can occur, depending on the angle 6 (see [3D]): 

if 9 is not of the form ^ for any integer k > 0, then 5 is reduced to a single point, and T has 
countably many multiple points, 

- if 9 is of the form ^ where k > 0 is an integer, then S is a Cantor set, and the Hausdorff dimension 
of S, noted dim# S, is dl ™ H 1 . 

2.1.3 The self-similar measure p. 

Recall the classical result on self-similar measures, see [IT, [21] and [25] page 26. 

Theorem 1. There exists a unique Borel regular probability measure p on T such that for any Borel set 

Act, 

= + ( 2 ) 

The measure /r is called the self-similar measure defined in the self-similar triplet (T, /i,/ 2 ). 

Let LP {T), p £ [1, Too) be the space of the measurable functions v on T such that / r |u| p dp < oo, 

endowed with the norm ||u|| iP ( r ) = (/ r |u| p d/i) 1 ^. A Hilbert basis of L 2 ^{V) can be constructed e.g. 
with Haar wavelets. 

The space W s,p (r) for s £ (0,1) and p £ [l,oo) is defined as the space of functions v £ L P (T) such that 
H w „,p(r) < oo, where 


M~<o - (/ r / r T^ 

Endowed with the norm IMI^.p#) = IM| iP ( rl + Mw s .p(r)> tlm spaces W S ’ P (T) are Banach spaces. In 
the special case p = 2, the space W s ’ p (r) is a Hilbert space, and is noted H S (T). 







The domains f2j nt and fl e xt 


Remark 3. In the special case when 9 = 0 and r = r@ = 1/2, the set T is in fact a line segment. This 
geometry will he discussed in Section\5\ (see Figure 01. In this case, it can be proved that if s € (0,1), 
then an equivalent norm of the space VF s,p (r) is given by 

IMlL P §'”(r) : = [ \ v f d h + J2 rkP [ \ v ~ ( v ) 2 U(r)\ P dx (3) 

Jr o *eA k j 2 f°( r ) 

where 2f a (T) is the intersection with T of the segment obtained by expanding the line segment f a ( r) with 
a factor 2 around its center (see 12PI (£$/)■ As in the rest of the paper, if v is a measurable function in a 
measured space (X,m), the notation (v)x refers to the mean value m (x) fx v dm - 

2.2 The domains f2 int and fl ext 

Call P\ = (—1,0) and P -2 = (1,0) and T 0 the line segment T 0 = [PiP 2 ]. Let us assume that /2(A) and 
/2(A) have positive coordinates, i.e. that 

rcosd < Pi and rsind < p 2 . (4) 

Let us also assume that the open domain Y° inside the closed polygonal line joining the points P\, P 2 , 
/2(A), /2(A), /i(A), /i(A), A in this order must be convex and hexagonal, except if 9 = 0, in which 
case it is trapezoidal. With dj), this is true if and only if 

(Pi — 1) sin 9 + p 2 cos 9 > 0. (5) 

Under assumptions a and a, the domain Y° is contained in the half-plane x 2 > 0 and symmetric with 
respect to the vertical axis x-\ = 0. 

Call K° = Y°. It is possible to glue together K°, fi(K°) and f 2 (K°) and obtain a new polygonal domain. 
The assumptions a an d a imply that Y° fl fi(Y°) = 0 and Y° fl f 2 {Y°) = 0- 



Figure 1: The construction of the first cell Y° 

Let the open domain fl int (see Figure a be defined as follows: 

U int = Interior | (J f a (K°)j , (6) 

Vo-e.4 / 

with the notations of H2.1.1I 

For a given 9, with rj defined as above, the following assumption on (a , p) will be made: 

Assumption 1. For 0 ^ 9 < 7r/2, the parameters a and P satisfy a an< 4 a f° r r = r o, an< 4 are such that 

{ i. for all r, 0 < r < rj, the sets Y°, f a (Y°), 0 € A, are pairwise disjoint, 
ii. for all r, 0 < r < r* e , /i(fA t ) D / 2 (flint) = 0, 

Hi. for r = r%, /1 (flint) H / 2 (flint) ^ 0- 
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Remark 4. Assumption^ implies that if r = rg, then /i(fli nt ) H — 0- 

In the case 9 = 0, Assumption [T] is satisfied by any a > rg = 1/2 and (3 > 0. The following theorem, 
proved in [5], asserts that for all 9 £ (0,7t/2), there exists {a, (3) satisfying Assumption 1. 

Theorem 2. If 9 € (0,7r/2), then for every a > VgCOs9, there exists (3 > 0 such that for all (3(^ (3, (a, (3) 
satisfies Assumption Q] 

Let D be an open bounded domain with a Lipschitz boundary, containing the closure of flint- The 
exterior domain fl ext is defined by 

flext := Interior (D \ fl int ). (7) 

Remark 5. The assumption that flint D may be relaxed: in fact, it would be enough to assume that 

r <e £>. 

Displayed on Figure [2] are examples of the domains flint and fl ex t, for the parameters 9 = 7t/5 in the 
left-hand side and 9 = 7r/4 in the right-hand side. 



Figure 2: The ramified domain fl for 9 = n/5 (left) and 9 = 7t/4 (right) when r < rg, f3\ = 0.7, (32 = 4. 


2.3 The truncated domain Q” t and the prefractal interface 

For every integer n ^ 0, the truncated domain f![( lt is defined by 

fl” t = Interior | (J (J / CT (A" 0 )J, (8) 

with the notations of 1 )2.1.11 As above, the exterior domain associated to f2^ lt is 

^ext = Interior(D \ ). (9) 

Note that the set T" defined in 1 )2.1.11 is a part of the boundary of f![( lt . The sets T n , n (2 0, will be 
referred to as prefractal approximations of the fractal set T. 

3 Function spaces 

Hereafter, we consider a domain fl; n t as defined in with 9 in [0, n/2) and r ^ rg, and we assume that 
the parameters a, (3 are such that Assumption |Tj is satisfied. 
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Trace results 


We define W 1,p (fl), p = [l,oo], = tti nt or fl = fi ext , to be the space of functions in L p (tt) with first 

order partial derivatives in L p ( 12). 

The space W 1 ' p (fi) is a Banach space with the norm (||u||£ p(fJ) + \\-§^\\ P LP(n) + 11^11^(0)) ^, see for 

example [6]. Elementary calculus shows that ||w||vyi,p(n) := (lMlz,p(n) + ll^ 7u lli,p(n)) is an equivalent 
norm, with ||Vu||£ p(n) := f n |Vw| p and |Vu| = y/\£^\ 2 + \-§£\ 2 ■ 

In the special case p = 2, the space W 1,p (12) is a Hilbert space, and is noted 77 1 (H). 

The spaces W 1,p (12 int ) as well as elliptic boundary value problems in fij nt have been studied in [3], with, 
in particular Poincare inequalities and a Rellich compactness theorem. The same results in a similar but 
different geometry were proved by Berger [5] with other methods. 


3.1 Trace results 

3.1.1 The classical definition of traces 


We recall the classical definition of a trace operator on du> when ui is an open subset of K 2 (see for instance 

[25 p. 206). 


Definition 1. Consider an open set oj C R 2 . The function u € L] oc (u) ) can be strictly defined at x £ oj if 
the limit 


u(x) 


lim 


\B(x,r) (~l w| y B (x,r)n^ 


u(y) d y 


( 10 ) 


exists, where \B(x,r) nw| is the 2-dimensional Lebesgue measure of the set B(x,r ) flu;. In this case, x 
is said to be a Lebesgue point of u. 

The trace is defined to be the function given by u\g u (x) = u(x ) at every point x € du ; such that the 
limit u(x) exists. 


Remark 6. Recall that for any p > 1, a function which belongs to W 1,P (R") can be strictly defined except 
on a set with zero p-capacity, see for example m and J26l 


3.1.2 A trace theorem on T 

It has been shown in [2] (see Theorem 11) that every function in W 1,p (12i n t) can be strictly defined on 
T 'H 1 -almost everywhere, where TL 1 is the one-dimensional Hausdorff measure. Moreover, the following 
trace result holds. 

Theorem 3. (see w 

• Assume r < r* s . For all p e]l,oo], if u £ W 1,p (12i n t), then W|r € W 1 ~’ P (T), and there exists a 
constant C > 0 independent of u such that 

H U l r llw 1-2 ^’ p (r) ^ C 'll u ll^ 1 . P (n int )- 

• Assume r = rj, then 

o the previous result holds if 1 < p < 2 — dim# S, 

o if p ^ 2 — dim# S, then W 1,p (f2i n t)|r C W s,p (r) for all s < l(d —dim#S), and the embedding 
is continuous. Moreover, if s > ±(d — dim#S), then W 1,p (12i n t)| r (jL W s,p (r). 

Remark 7. The space of the traces on T of functions in bE 1,p (f2j n t), 1 < p < oo was characterized in w> 
whether r < r# or r = r* s , as the space JLip(l — ^-^-,p,p;0;T), which was first introduced in [23f. Of 

course, if r < r* s , then JLip(l — ^p,p,p;0;T) coincides with W 1 ~’ p (r). An easy consequence of this 
characterization is that the space of the traces on T of functions in W 1,p (f2 int ) is relatively compact in 

L p ( r). 
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Proposition [I] below will be useful in the proofs of the main theorems of this paper. 
Proposition 1 . For every u £ H 1 { fli nt ), 





Proof. The present proof relies on Proposition 1 in [3], which states that for any u £ H 1 ( fli nt ), the 
sequence of piecewise constant functions (u n )neN defined on T: 


u n — ( M )r^ 

cr&An 


where (u) r „ = J r „ Ujr^(z) d z, is such that limn._j.oo \\u n — U|r||,L 2 (r) = 0. Note also that f r u^d/.i = 
Ea^A jr^j ( u )r tr - Hence, in order to prove Proposition [TJ it is enough to prove that 


pq f rn u |r- 2 da; - f r u 2 n dp. 


0, or in an equivalent manner, that S n 


0, where 


Sn = pij U M Z ) d " - l r<T ! <«>£-) = jj^j [ u \r°(z) - («) 


dz. 


From a standard trace result on T° and appropriate rescalings, we know that for a positive constant 
independent of n, o £ A n and u £ 7J 1 (17i nt ), /r- [u| r <r( 2 ) - (u) r „] dz < C\T a \ /^ ( n lnt) |Vit| 2 . Hence, 

S n < CEaeA„ pf //„(n lnt) l Vu l 2 = % I fA n int) l Vu l 2 ’ which im P lies that s n °' D 

We also recall the following refined trace inequality, we refer to [3] for the proof. 

Theorem 4. [see EP. Th. 11] Assume that r ^ 1/2, then for all real number n £ (2r 2 ,1), there exists a 
constant C such that for all v £ H 1 ( fh nt ) with z>|po = 0, 

hrllljp, « C •£ «” E ll^ t, llL 2 (/ r (y°))- (11) 

ra^O r£An 


3.2 Extension results 

3.2.1 The subcritical case r < r$ 

As seen in § [1] it was proved in 0] that if r < rg , then flj nt is an (e, d)-domain (see [22] for a definition), or 
in an equivalent manner, a quasi-disk (see [32]). Hence, the extension result of Jones and Yodop’janov et 
al. applies, and fli nt is a Sobolev extension domain (see [22), be. Hi nt has the VF 1,p -extension property 
for every p £ [1, oo]: there exists a continuous linear operator £ from kF 1,p (Hi nt ) to W 1,P (R 2 ) such that 

S(u)\n int =u, \/u £ W 1,p (n int ). (12) 

Similarly, the set H ext = D \ fij nt is an (e, d)-domain, and thus a Sobolev extension domain. 

3.2.2 The critical case r = r^ 

When r = r%, it is easily seen that fh nt is not an (e, d)-domain, and the extension results of Jones and 
Vodop’janov et al. do not apply. In fact, if p £ (1, oo), it is easy to construct a function u £ W l,p {£l m t) 
such that u = 1 in /i(fl; n t) and u = —1 in //(Hint). If p > 2, u cannot be extended to a function 
belonging to JF 1,P (R 2 ) because the existence of such an extension would contradict the Sobolev imbedding 
of W 1 ’ P (R 2 ) in C(R 2 ). 

In the case when r = rj, the situation depends in fact on the Hausdorff dimension of the set E = 
/l(r) n/ 2 (r). The following extension theorem holds. 
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Theorem 5. see TlB L [2| / Set p* = 2 — dim# 2 (recall that cl is defined in (jT|) ). 

1. If p £ (l,p*), then flint has the W 1,p -extension property. 

2. If p > p*, then fl; n t does not have the W 1,p -extension property. 

Point 1 in Theorem [5] was obtained in m • Point 2 is a consequence of [I] and [2]: by Theorem [3] if 
p > p*, then hP 1,p (riint)| r TT 1 ”~’ p (r) = W 1,p (M 2 )| r . This is in contradiction with the existence of 
a continuous extension operator from W 1,p (fli n t ) to 1V 1,P (R 2 ) (see [1^ for the proof that the notions of 
traces coincide /z-almost everywhere on T). 

Remark 8. ^4s it was seen in $2.1.21 only two situations can occur, depending on the geometry ofH lnt : 

- if 9 is not of the form w/(2k) for any integer k, then dim# 2 = 0 and p * = 2, 

- if 9 is of the form tt/( 2k) for an integer k, then dim# 2 = (dim# r)/2, and p* = 2 — (dim# T)/2. 

Remark 9. The special case p = p* is not dealt with in Theorem [31 The latter is of particular importance 
in case 1 of Remark 0 above, since the case p = p* = 2 corresponds to the question of the H 1 -extension 
property. 

In fact, it was proved by Koskela in wa that if a domain in R™ has the W 1,n -extension property, then it 
must have the W 1,p -extension property for every p ^ n. Hence, a consequence of Theorem [5] is that fli nt 
cannot have the W 1,p -extension property when p = 2. In particular, the domains that we will consider in 
Section^ fail to satisfy the H 1 -extension property. 

To the best of our knowledge, the question of the extension property for p = p* in case 2 of Remark [21 
seems to be open. 


4 The transmission problem in the case r < 

4.1 The transmission problem with fractal interface 


The 


transmission problem can be 




formally stated as 


-Au = f 

in fli n t 

[u] =0 

on r, 

[d n u\ = au 

on r, 

int = fin^leyi t = 0 

on E, 

fin^leyit — 0? ^int — ^0 

on r°, 

fin^leyit = 0 

on dD 


(13) 


where a > 0, E = <9fli nt \ (TUT 0 ), [u] (resp. [9„u]) denotes the jump of u (resp. of the normal derivative 
of u) across T, f £ L 2 (D) and uq £ H 1//2 (r°). We also use the notations u; n t •'= W|n int and u ext := u |n ext - 
The transmission condition [d n u\ = au on T has no real meaning, since the normal is not defined on T. 
The rigorous meaning of m is the following variational formulation: 


find u £ V such that for all v £ Vo, 


a(u, v) 


where V is the affine space defined by 


fv dx, 


id 


(P) 


V" — {ti £ T“(H), Ui n t £ H (fli n t); ^ext £ H (f^ext); ^int|r° — ^0; ^int|T — ^ext|r}- (14) 


Recall that v £ >->■ U| r £ H ll/ 2 (r) and v £ H l (i l ex t) '-t U| r £ iJ 1 / 2 (T) are continuous maps, 

hence V is closed. Note that if u £ H 1 (D) and = u o then u £ V. 

The vector space Vo is defined as V, except that the condition Ui„t|ro = uq is replaced by Ui n t|r° = 0. 
Finally, 


l ci in 


f 

Jn 


Vu int • Vc int dx + V u ex t ■ Vu ex t dx + a U|r^|r dp. 


a(u, v ) 


(15) 
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Remark 10. The traces in the condition Ui n t|p = u ex t|r (El) are meant in the sense of Definition Q] 
The above definition of the space V is suitable when r < r@ and in the special case when r = r@ and 6 = 0 
discussed in Section 0 In the other cases, the transmission condition has to be considerably changed, see 
Remark m 

Note that the space Vo, equipped with the norm a(w,tt) 1 / 2 is a Hilbert space. From the Lax-Milgram 
theorem, we see that for every function / given in L 2 (D), there exists a unique weak solution u G V to 
(P). Moreover, u minimizes the functional 

v G V a(v, v) — 2 / fv da;. (16) 

JD 


4.2 The transmission problem with prefractal interface 


For any positive integer n, let us consider the similar transmission problem in which the interior domain 
has been truncated by stopping the construction at step n. This class of problems is much more standard 
since the interface T n consists of 2 n pairwise disjoint line segments. The boundary value problem reads: 


-Ait = / 

in 

^nt 1 

M = o 

on 

pn 

,e, 

II 

3 s 

s 

on 

pn 

dnK lit = S «< xt = 0 

on 

E", 

9nK xt = 0, «£ t = Uo 

on 

r°, 

dnU^t = 0 

on 

dD, 


un: 


(17) 


where £” = <912^,. \ (r"ur°), and [it] (resp. [d n u]) denotes the jump of u (resp. of the normal derivative 
of u) across T". We also use the notations u-^ = u^n and u" xt = 

The variational formulation of problem m can be stated as follows: 

find u G V n such that for all v G Vq 1 , 

a n (u,v ) = / fvdx, 

J D 

where V n is the affine space defined by 

V n = {u£ L 2 (D), < t G < xt G < t|r o = u 0 , < nt|r „ = < xt|r „}. 


(Pn) 


(18) 


Remark 11. Let G n be the closure of the set 90^ \ (r° UF 1 ), which is a finite union of polygonal lines. 
It is easy to see that V n is the set of the functions in H 1 (D \ G n ) such that Wint|r° = u o- 
Similarly, if we define G as 9Sli nt \ (r° LIT), we observe that G is not closed, since its closure contains 
r. Observe that, in general, the functions u G H 1 (D \ G) do not satisfy Ui n t|r = M ex t|r, so V cannot be 
identified with the set of the functions u G H 1 (D \ G) satisfying the Dirichlet boundary condition on T . 
On the other hand, since D\G is not an open set, dealing with H l (D \ G) is not very straightforward. 

Here also, VJJ 1 is defined as V n , except that the condition = «o is replaced by 'u[j lt , r o = 0, and 

a n is defined by 


a n (u,v) = 


V< nt -V< nt dx 


Vu P , 


VCt dx + 


a 

jf" 


M|rnV|pn dx. 


(19) 


The space Vq , equipped with the norm a„(u, it) 1 / 2 is a Hilbert space. We also remark that V n C V 
with a continuous imbedding. Again, the Lax-Milgram theorem implies that for every function / given 
in L 2 (D ), there exists a unique weak solution u n G V n to that problem, and u n minimizes the functional 


v£V n i->- a, 


i{v,v) -2 f fv dx. 
J D 


( 20 ) 
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Lemma 1. The sequence u n is bounded in V. 

Proof. Let u o be a function in H 1 (D) such that Mo|r° = wen and such that uq is supported in a compact 
set which does not intersect the sets T", V?i >1. It is clear that uq G V and that uq € V n for all n> 1. 
Let us define Co = f D (|Vito| 2 — 2/fto) dx. Thus, for all n > 1, 

a n (u n ,u n )- 2 / fu n da; < a n (u 0 ,u 0 ) -2 / fu 0 da; = C 0 , (21) 

J D 

because u 0 |r™ = 0. On the other hand, since V n C V, 

a n (u n ,u n ) ^ / |Vu n | 2 dx + / |Vu n | 2 dx = / |Vit n | 2 dx + / |Vu n | 2 dx. (22) 

J J Qg xt ^int J ^ext 

We shall also use the following Poincare inequality: there exists a constant C > 0 such that, for all v G V, 
Mh(D) < c (lhr°||| 2 (r°) + f |Vu| 2 d x+f |Vu| 2 dxV (23) 

\ J £^int J ^ext / 

From (l2ll) . (l22l) and (l23l) . we deduce that 

1 

f |Vu„| 2 dx+ f |Vu n | 2 dx-2\/C||/|| i 2 (I j ) 

J ^int J O ex t 


ll^oir 0 IIz, 2 (r°) + j |Vu n | 2 dx + j |Vi 

J r^int j ^ext 


dx 


^ Co 


which implies that the quantity f n |Vm„,| 2 dx + f n |Vw ra | 2 dx is bounded by a constant independent 
of n. Using (l23l) again, this implies that ||uri||z, 2 (D) is also bounded by a constant independent of n. 
Combining the previous two observations, we obtain that the sequence u n is bounded in V. □ 


4.3 M-convergence of the energy forms in the case r < r* e 

We start by extending the definition of the forms a and a n to the whole space L 2 (D) by setting 

a{u, u) = oo if u £ L 2 (D) \ V, 
a n (u, u) = oo if u € L 2 (D) \ V n . 


(24) 

(25) 


We are interested in proving the convergence of the forms a n to a in the following sense, introduced by 
Mosco (see [55]). 

Definition 2. A sequence of forms (a„) is said to M-converge to a form a in L 2 {D) if 
i. for every sequence {u n ) n weakly converging to a function u in L 2 (D), 

lim a n (u n ,u n ) ^ a(u,u) as n —> oo, (26) 

ii. for every u € L 2 (D), there exists a sequence {u n ) n converging strongly in L 2 {D) such that 

lim a n (u n ,u n ) ^ a(u,u) as n —> oo. (27) 

Theorem 6. Assume that r < Vg, then the energy forms a n M-converge in L 2 (D) to the form a. 

Remark 12. The M-convegence of forms differs from the Y-convergence only in that the sequence ( u n ) 
in point [I] of Definition^ is assumed to converge weakly instead of strongly. In the following , only the 
Y-convergence of the energy forms a n will be needed. 

Proof. We will prove separately points [T] and [IT] in Definition [2] 

Proof of point\Y\ Suppose that (u n ) weakly converges to u in L 2 {D). Without loss of generality, one can 
suppose lim a n (u n , u n ) is finite. We may further assume that there exists a subsequence, still called (u n ), 
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such that a n (u n , u n ) converges to some real number as n —> oo; as a consequence, there exists a constant 
c independent of n such that 

a n (u n ,u n ) ^ c. (28) 

In particular, for all n, u n £ V n 1 which implies that u n £ V. Then, (1661) implies that (un|n int ) is 
bounded in _ff 1 (rii n t), and (wn|n ext ) is bounded in i/ 1 (f2 ext ). Therefore, there is a subsequence that we 
still note ( u n ) such that ('(tn|n int ) converges weakly in H 1 (f2i n t), and strongly in L 2 (fli nt ). Since u n —*■ u 
in L 2 (D), we see that Wn|n int —*■ Mint in t ). Similarly, up to a further extraction of a subsequence, 

u ™|fiext -s ‘ Uext in ^(^ext)- Consequently, 


lim 


( f |Vu ra | 2 dx+ f |Vu n | 2 dxj ^ f |Vu| 2 dx+ j |Vit| 2 dx. 

V^Sit / O j nt J&ex t 


We will now prove that p4r f pn u 2 da; —> f r u 2 d fi as n —> oo, which will yield point [T] 

The following inequality was proved in [3]: for every v £ H 1 ( flint), ||u — (^)r° ll^ 2 (p) ^ t )- 

This implies that 

|(u) r o| < IMUj(r) + C'l|Vn|| i 2( n . nt >. 

Similarly, ||u— (u) r o || i2 ( r o) < ^11 Vu|| i2 , n . , implies that 

yj^lMU a (r°) < Kt>) r o| + C||Vt;|| La(nint) . 

for some constant independent of v that we still note C. Combining these two inequalities, we obtain 
that 

: IMU 2 (r°) < IMUj(r) + c 'l|Vf|| L 2( n . nt ). 




Hence, for every a £ A„ 




IMk 2 (r-) < ||uo/ CT || i 2 (r ) + C'||V(uo/ CT )|| i2(nint) 

= 2 2 ||u|| L 2 (/iT(r)) +C||Vu|| i2(/(r(nint)) . 


This yields that 
1 


< 2 n |Mlz, 2 (/„(r)) + 2C2 2 IK’IUj(/ CT (r))l|Vu|| i2(/(7(nint)) +C ||V«|| L2(/<T(fllnt)) . 

Therefore 


\ T r, 


i dx = y — 

1 2 n ^ |T CT 

O-eAn 1 


■Ur, 2 da; 




C 2 ■ 9 

|| 2(r )+ 2C2“ 2 ||w„|U2 ( /„ (r)) ||Vu„|| i2(/t7(n . nt)) + — l|Vu„||22 (/iT(f2int)) 

O-eAn O-eAn 




c 2 


\u n \\ 2 L 2 tT) + 2C2 2 ||u n || L 2 (r) ^2 ll Vu «lli 2 (/„(n int )) + ll Vu "lli 2 (.Mn iri t))' 




O-eAn 


Since u„|Q. nt is a bounded sequence in H 1 ( flint), there exists a constant M such that 


jp^T J u n 2 da; ^ IIIIl 2 (r) + 2CM2 2 ||un||Lj(r) + 


C 2 M 2 


Moreover, since u n weakly converges to u in H 1 (Hi nt ), then up to the extraction of a subsequence, u n |p 
strongly converges to U|p in L 2 (T), from Remark^ Hence, we obtain that 
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Similarly, the following inequality holds for every v G H 1 (flint): 

IMIl=(T) ^ ^jj^j-IMlL 2 (r°) + C H^' l ’llL 2 (n in t) (29) 

for some constant C independent of v. As above, we deduce that 



and we obtain the desired result. 

Proof of point \nf\ Take u G L 2 (D). By (l24l) . we may assume that u G V. We must construct (u n ) 
converging strongly in L 2 (D) such that E71) holds. Note that the choice u n = u, cannot be made, since 
u “fL V n in general. 

Take 6 > 0 and consider a neighborhood u C D of flint such that flint <s to and sup xGul d(x, flint) < 8, 
where d(x, flint) = inf ye n lnt |x— y\- We introduce the notations to a = / CT (w) for a £ A and ui n = Ucre^t W<T 
for all integer n. 

For every n, introduce the cut-off function \ n in D defined by 

Xn(x) = (1 - 5r~ n d(x,u] n )) + , (30) 

where a + stands for the positive part of a real number a. Hence, \n = 1 in to n and Xn = 0 outside 
Co n := {x G D , d{x,uo n ) < Sr n }. Note that if we set Co := {x G D , d(x,co ) < <5} and Co a := f a {Cj) = {a: G 
D, d(x,oj a ) < <5r™} for a G A n , then Co n = LLeAi ‘ 

We can assume that S is small enough so that Co a D Co T = 0 when cr, r G A„ and cr ^ t, since / (T (flint) FI 
fr (Flint) = 0- 

We now define a sequence of functions u n by 


u n = {\- Xn)u + Xn£{u int), 


(31) 


where £ is the extension operator introduced in m and as above, u; n t = U|n int - Obviously, u n belongs to 
the space V n and the sequence (u n ) strongly converges to u in L 2 (D). We will prove that lima„(u n , u n ) = 
a(u, u) as n —> 00 . We start by showing that 



Hence, it is enough to show that f- , t |V(u — u n )\ 2 dx —> 0 as n —> 00 . Note that 




\\7(Xn(£{uint) 


-u))\ 2 dx ^ 2 (if + lf) 


(33) 


where If = I V(£(u int ) - u )| 2 dx and I 2 = |(Vxn)(£(« int) - u )| 2 dx. 

First observe that if — > 0 as n —> 00 since V(£(iti n t) — u) G L 2 (fl ex t)- We are left with dealing with If . 
One has 

if < cr -2 " f \£{u int ) - u\ 2 dx 

J LOn\Qint 

< cr~ 2n ^2 ( l^(uint) - u\ 2 dx. 

O'^'A.n 
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where c > 0 is a constant independent of n. Introduce the set 6j = f\ 1 (/i(w) \ flint)- We have the 
following Poincare inequality: there exists a constant C such that for every v £ H l (w) such that U|p = 0, 



d x^C |Vv| 2 dx. 


(34) 


Observe that if S is small enough, then cD CT \ f l int = f a ( cD) for every a £ A. Therefore, applying a rescaled 
version of © to the function £’(ui nt ) — u, we obtain that there is a constant c' > 0 indenpendant of n 
such that 




d ^2, f |V(£(u int ) - u)\ 2 da: 

CT^A-n ^ \^int 

c' f |V(f (u int ) — u)| 2 dx 

^ \^int 


since the sets uP, a £ A n are pairwise disjoint. We deduce that 1% —> 0 as n —> oo, since £{u- m t) — u £ 
H 1 ( fl ex t), which yields (15^1) . 

We will now prove that 


|r n 


J u n 2 dx — J 


dfj, 


(35) 


as n —> oo, which will conclude the proof of point \n^ Observe that for every integer n, £(uint)|r" = u, 


which implies that u n | r » = u. We are left with proving that j^n- f r „ u 2 da; 
which holds by Proposition [T| □ 


f r u 2 d/i as n —> oo, 


A standard consequence of the Mosco-convergence of the energy forms proved in Theorem [G] is the 
convergence of the solutions of the problems ( |P„| ) to the solution of problem |P]), in L 2 (D) and in V 
(recall that u n is bounded in V from Lemma [I]). 

Theorem 7. Take f £ L 2 (D), and note u n (resp. u) the solution of problem fl_P n | ) (resp. (fP)) ). The 
sequence (u n ) converges to u in the space V. 


5 A particular geometry with r = 

As seen before, the proof of Theorem [5] is based on the extension result of i J3.2.1l In the case r = rj, the 
-/^-extension property is no longer true for the domain f2i nt (see Remark [9]). In what follows, we focus 
on the special case when 9 = 0. We will see that in this case, the transmission condition imposed on T 
yields an extension result (see Theorem H) which is the main ingredient for proving the M- convergence 
of the energy forms. 

In the case 0 = 0, it can be seen that rf = and the ramified domain described in m is as in Figure 
[3j In this particular case, the set /i(r) D / 2 (b) is reduced to a single point that we call A. Observe that 
the self-similar part T of the boundary is a line segment, and the self-similar measure p, associated with 
r is the normalized one-dimensional Hausdorff measure. 

Since r = rg, the domain O ex t has infinitely many connected components. Call U the outer connected 
component of O ext , ie the only connected component which has a nonempty intersection with dD (see 
Figure [3]). Observe that T is a subset of dU, and that the intersection of T with the boundary of every 
other connected component of O ex t is reduced to a single point. Apart from U, each connected component 
of flext is a triangle whose top vertex is at a dyadic point of T. The largest triangle is named P, see 
Figure^ and all the other triangles are the images of T by f ai a £ A n , n > 1. 

We consider the transmission problem 

—A u + (3u = f 

^int|r = ^ext|r 
[d n u\ = au 

int = ^n^ex t = 0 
^n^ext — 0? ^int — ^0 
k ^n^ext = 0 


in flint u a 


int ^ ^ £ ext 5 


on r, 

on dfl int \ r, 
on r°, 
on dD , 


( 36 ) 
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where a and f3 are positive numbers. The trace rt ex t|r i n the transmission condition is meant as the trace 
of the function u\jj on the set T. 

Remark 13. The reason for considering the operator —A u + (3u with fj > 0 instead of —Au as in the 
former case is that, in the present case, f2 ext is an infinite union of disjoint connected sets: !7 ex t = 
uu{j aeA faiT). Therefore, 136 1) involves Neumann problems in T and in fa(T), a £ A n , n > 1, which 
are not well posed if (3 = 0 and the average of f in these sets is not zero. It would also be possible to 
consider the case (3 = 0 under additional assumptions, on the support of f for example, but this would 
imply further technical details, because the solutions of Neumann problems in the holes would then be 
defined up to the addition of constants. 

Remark 14. When 9 > 0 and r = r*, the situation is quite different: T is not entirely contained in the 
boundary of any connected component of fl ex t- It can be shown that there exists 6 £ (l,d) such that the 
intersections of T with the boundary of the connected components of fl ext have Hausdorff dimension 8. 
These sets are called the canopies of the domain f2j nt . In this case, the transmission condition has to be 
described more carefully. This is the topic of a work in progress. 

The meaning of (Idfil) is the variational formulation ([21) where V is defined by (TT4l) and 

a(u, v) = / (Vuint-Vuint+zStiintUint) dx + / ('Vu ext • Vu ext + 0u ext v ext ) dx + a / U|rU|r dp. (37) 

J Qint I Q ext Jr 

In order to set the transmission problems in the geometries with prefractal interfaces, we first need to 
define some trapezoidal subsets of the triangular holes as follows: let H be the height of the triangle T. 
Choosing the coordinates in such a way that T° is a segment of the line {X 2 = 0}, we see that T is a 
segment of the line {X 2 = 2 H}. Then, we can also define T n and T” by T n = Tfl {(2 — 3/2 n+1 )H < 
x 2 <{2- 2 ~ n )H} and ff = f a {T) (~l {(2 - 3/2 n+1 )H < x 2 < (2 - 2 ~ n )H} for cr e A rn and m < n. 
Finally we define 

n —1 

^xt = U U <= ^ext n 

m —0 a^Am 

The transmission problem with interface T” is then 


2H — 


2 n+1 


H < x 2 < 2H - 2 ~ n H 


(38) 


where 


—div (v n Vu) + f3u = f 

in ^Pnt 

II 

0 

on r n , 

3 

II 

T 

2? 

on T 71 , 

9nK Zt = 9 n < x t = 0 

on E n , 

<9n< xt = 0, M/nt = U 0 

on T°, 

e 

c+- 

II 

0 

on dD 

«/n = 2- 2 "l S n t +l 



(39) 


(40) 


Note that v n = 1 in 0/: lt and that the Lebesgue measure of the set where v n = 2 2n vanishes as n —> 00 . 
Hence v n tends to 1 almost everywhere in D. The variational formulation of (l39l) is with V n defined 
in m , and a n defined as follows: 


a n (u,v) 


/ (V< t • + /?<„ t <t) dx + / (^„VUe Xt • V< xt + /3< xt i4 t ) da; + / U| r nW| r n dx. 

/nPn. l r I • /r " 

(41) 


Remark 15. The reason for modifying the partial differential equation in \39\ ) by taking — div(i'nVu) 
instead of —Au in m is that for a function u £ V, rt|o int is completely independent from U|/ ct (t)- This 
explains why the construction of a sequence of functions (u n ) such that u n £ V n , u n —»• u in L 2 (D) and 
a n (u n , u n ) a(u,u), is difficult without modifying the coefficients of the partial differential equation near 
the top of the triangles T and fa(T) in order to cope with the possibly strong gradients of u n . Although 
we have not tried it, it may be possible to choose a parameter larger than 2~ 2n in the definition of v n . 
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D U 


Figure 3: An exemple of the domains Qj n t and fi ext in the case 0 = 0 


The main result of this paragraph is the following theorem 

Theorem 8. Assume that 0 = 0 and r = r] = [ . Then the energy forms a n defined in EZP M- converge 
in L 2 (D) to the form a defined in 

Since dU is Lipschitz-continuous, a standard trace result yields that for every u £ V, w e xt|r £ H l / 2 ( r). 
Hence, the transmission condition in (1361) implies that 

Mint I r £ H 1/2 (V). (42) 

Note that (TE1) is not only a consequence of the fact that iti nt £ -H 1 (^int), because the latter property 
only implies that Wi n t|r in H S (T) for all s < \ (see Theorem[3|. 

For proving Theorem [5J we need the following extension result, which is not available in the literature: 

Theorem 9. There exist a linear extension operator F from {y £ H 1 (fl j nt ), t>| r £ H 1/l2 (T)} to H 1 (D) 
and a constant C > 0 such that for every v £ H 1 (fl i nt ) with z>|r £ H 1 / 2 ( F)}, 

IHI.fr 1 (z>) ^ + lk| r llif i / 2 (r)) ■ (43) 

Theorems [Till andfTTlbelow will play an important role in the proof of Theorem[9] We start by recalling 
an extension result for multiple cones from [7]. 

Theorem 10. [see B Th. 5.1] Call C the double cone in R 2 defined by |xi| < \x 2 \. Write p(x) = ||x|| 
(the notation || . || stands for the euclidean norm), and, for every v £ H l (C), introduce the antiradial 
part v a of v in the cone C, defined by 


v a (x) = v{x)-(v) Sp(x) , 

where for any R> 0, Sr = {x £ C, ||x|| = R}, and ( v) Sr is the mean value of v along the arc Sr. 
There exists a linear extension operator 

A : jv e H\C), ^ € L 2 (C) | H\R 2 ), (44) 

such that for every v £ {u £ i7 1 (C), v a /p £ L 2 (C)}, 


II At 


\m( 


2) ^ c Ipllffi(c) 


LHC); 


where c > 0 is a constant independent of v. 
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Remark 16. The construction in E is such that if v is radial (resp. constant) in C C1 5(0,1?), then Av 
is radial (resp. constant) in 5(0,??). 

As mentioned in [7], TheoremQIJlcan be immediately extended in R n to the case of a union of two half¬ 
cones sharing the same vertex, separated by a hyperplane passing through the vertex and not containing 
any direction of the boundaries. 

Theorem 11 : Peetre-Tartar. [see S3 El/ Let E, E\, 5 2 , F he Banach spaces, and let Ai, i = 1,2, be 
continuous linear operators from E to Ei, and suppose A\ is compact. Further assume that there exists 
a constant cq > 0 such that for any v £ E, 

IMI e < c 0 (||Aiu|| Ei + p 2 u|| B2 ). (45) 

If L is a coninuous linear operator from E to F such that A|kerA 2 — 0, then there exists a constant C\ > 0 
such that for any v £ E, 

II-MIf < ci||A 2 u|| S2 . (46) 

Notations We start by introducing notations for the proof of Theorem [9] 

We first introduce a domain C which is the union of two truncated half-cones included in flint, whose 
common vertex is the point A. Recall that T is the main hole of the domain flint- Call ipo £ (0, \) the 
upper half-angle of the triangular domain T (see Figured]), and take > ipo- Call C the half-cone whose 
boundary is made of the two half-lines through A with respective angles ipo and ipi with the vertical axis 
(see Figured]). Call C 2 = CD (flint \T°). We can assume that <pi > ipo is small enough so that C 2 C flint, 
in other words C 2 does not intersect any of the holes of fl; n t- We define C\ to be the symmetric of C 2 
with respect to the vertical axis x 2 = 0, and we write C — C\ U C 2 . 

We also introduce the sets Y k ’ 1 := /1 o f 2 k ~ 1 (Y°) and Y k ’ 2 := / 2 o f 1 k ~ 1 (Y°) for every k > 1 (see Figure 
|3|, and we write Y k := Y k ' x U Y k ' 2 . We also note 7 := /i 2 (r) U / 2 2 (r) (see Figured]). 



Pi O P 2 


Figure 4: The region C = C\ U C 2 and the domains Y k ’ 1 , k,i= 1, 2. 

Call fl 1,1 := /1 (flint \ y*2(f^int)) and ^l 1 ’ 2 : = /2(flint \ /1 (flint)) (see Figured]). We introduce the sets fl fc,i 
defined by fl fc ’* := <? fc_ 1 (fl l ), k ^ 1, i = 1,2, where g is the homothety centered at A with ratio 1/2, see 
Figured] 

For every integer k ^ 1, we introduce £ fc := {r £ A, f r (Y°) C fl fe,:L U fl fe ’ 2 }, with the notations of 1 )2.1.11 
Take k £ (1/2,1). We introduce the space G = (v £ L) oc (fl 1 ’ 1 U fl 1 ’ 2 ), |M| G < 00 }, where 


\\vf G = 


E- 

1 


E 

T€E 1 r\A„ 


lul dx. 


lfr(.YO) 


(47) 
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7 = /i (r) u/ 2 (r) 


Figure 5: The domains k,i = 1,2 (in light grey: ST 1 ’ 1 and f2 1 ’ 2 , in dark grey: fl 2,1 and fl 2,2 ) and 
the sets /) 2 (r), i = 1,2. 


Endowed with the norm || • || G , the space G is a Hilbert space. 

We also introduce the space H = {r; £ L 2 (H 1,1 U fl 1,2 ), Vv £ G}, which is a Hilbert space with the norm 

(l|Vr'|| G + ||U|| 2 2(Q1,1 U Q1,2)^ • 

Moreover, from Theorem U we see that v £ H H > U | 7 is a continuous operator from H to L 2 ( 7 ). Arguing 
by contradiction, we can show that 

is an equivalent norm on H. 

We first state and prove two lemmas which will be useful in the proof of Theorem [9] 

Lemma 2. There exists a constant c > 0 such that for every v £ H, 

{v) Yl ,!\ 2 + \v(x)~ {v) Yl , 2 \ 2 ^j dx < c(^J\v~ (u } 7 | 2 dn+ ||Vv|| G ^ . (48) 

Proof. We introduce the Hilbert space E = {(v,w) £ H x T^(y), U | 7 = w}, endowed with the norm 

II • \\ E g iven b y IIK w )IIe = II Vw IIg + IMIljo)* 

We now introduce the operators 


M ■■ {v,w) ££4 (u) 7 , 

A 2 : (v,w) £ E >->■ (Vv,w- (w) 7 ), 

L : (v,w) £ E (v - (v) yl:1 ,v - (v) yl:2 ). 

With the notations of Theorem fill E\ = R, E 2 = G x L 2 ( 7 ) and F = L 2 (F 1 ) 2 . It is easily seen that Ai, 
A 2 and L are continuous linear operators, and A\ is compact. Moreover, © is clearly satisfied. 
Observe that (v,w) £ E lies in ker A 2 if and only if v is constant in fi 1 ’ 1 and in f2 1,2 . Hence, it is obvious 
that L iker .42 = 0- From this, we deduce by Theorem [TT] that there exists a constant c > 0 such that 
\\L(v,w)\\ F < c||A 2 (u, w)|| B2 for all (v,w) £ E, which yields (1151) . □ 

Lemma 3. Assume that v £ H 1 ^; nt ) and U|r £ iF 1 ^ 2 (T) ; then ^ £ L 2 (C), and 


L 


dx 


< 


c 



M-J 2 <fa + 


E KP ~ k+1 

k^lp^k 


E 

T£?. k C\A„ 


IfrlYO) 
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for some constant c > 0 independent of v, where 7 * = g k 1 ( 7 ) (recall that g is the homothety centered at 
A, with ratio 1/2,). 

Proof. We first observe that, by self-similarity, there exists a constant Ci > 0 such that for all x G Y k , 
p(x) ^ ffr. Therefore, 

2 

f Hi dx < Cl V2 2k f \v a \ 2 dx : 

since C C U fc >i Y k by construction. Hence, there is a constant c 2 > 0 such that [ 

" Jc 

h = 5 Z 22fc [ \v(x)-{v) Yk \ 2 dx, 
k >1 jYk 

J 2 = J2 22k f “ MsrfJ 2 da; - 


where 


f 

V a 

lc 

P 


dx < c 2 (Ii+h), 

(49) 

(50) 


We start by dealing with I\. We note that 


h < 


f (\v-(v) Yk .A 2 + \v^(vM 2 ) 

2 k> 1 Jyk V 7 


dx. 


(51) 


For every k ^ 1, we can apply Lemma [2] to the function v o g k 1 . Since = g k 1 (y), we obtain 


h < c (^ 2 fc f \v- (u) 7 J 2 d M + ^ 

rGSifl. 


r\A m J 9 k ~ 1 (U(Y 0 )) 


\\7v\ dx 


— c 


E 2 7 I»-(»>JW + EE '= P_ '‘ +1 5 : 


ife>l 




|Vw | 2 dx 


(52) 


for some constant c > 0 independent of v, since {g k 1 o f T , r G E 1 n A m } = {/ T , r G Y, k D A m +k-i}- 
Let us now deal with J 2 . For every R > 0 and i = 1,2, call S l R = SrC I C, and (x) S i the mean value of 


on the set Sjj. We note that 


2 = 1,2 w y 7 


(53) 


Take z G {1, 2} and Xo G fi(Y°), and, for every integer k ^ 1, pk = p(x o)/2 fe 1 . 
We observe that Yk 2 2fc /yfc |(u)y*,i - (^)gi | 2 < 2(Ji + J 2 ), where 

p(®) 


Jl 

= ^ 2 2 fc + 1 |F^| {{v) Yk ,-{v) sip )\ 

(54) 


k^l 


J2 

= J2 22k [ l^s* — <«)s* I 2 da:- 

z ' Ivk Pk P(n) 

Z-^.1 J r 

(55) 


Let us first examine J\. The following Poincare inequality holds in Y 1,1 = fi(Y °): for every x G iL 1 (F 1 ’*), 

f |x(x) — (x) s | 2 dx < M f |Vx | 2 dx, (56) 

Jy !.* Pi 


where the constant M > 0 is independent of x. 
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Observe that for every integer k > 1, Y k ’* = g k 1 (F 1 ’ 1 ), and S l pk = g k 1 («S'p 1 ). Then 

1 ' 


\{v) Y k,i (v) S i 


/ (v(x) - (v) S i ) dx 

J Y k,i p* 








| Y k 

w^\l Y j v{x) - {v)si J dx 

2^-l\Y k ^\ Xm |V ° 9k ~ 1(x) ~ {V ° 9k ~ 1)s « |2 ^ 

2 2(fe-]0|yfc,i| J Ylii l V(t, °g fc ~* 1)|2 dx 


M 


2 2(k-l)\ Y k,i\ J yk i 

where we used (l56l) . From this and (l54l) . we deduce that 


|Vx | 2 dx, 


Ji < 8MV / |Vx| 2 . 

I Y k ’ i 


(57) 


k> 1 1 


To deal with J 2 , we use polar coordinates (p, ip) centered at A. Introduce the positive constants Rq,Ri 
such that for all ( p , tp) G Y 1 ’ 1 , R 0 < p ^ i?i. Observe that if x G 


Ms*,., “Ms* I 


951 


pi 


1 

- / MpOe), <p) - v(pfe, <p)) dtp 

<P0 y (^0 

— P P^^W) dsdtp 

— ^0 J Pk dp 


' VO J Pk 

rV 1 CPO) 






rvi r 

J<PO J Pk 


(¥>1 - T’o) 2 

/*<£l /* 

- ^0 Jip 0 jR 0 /2 k - 1 

lo s§ 


W Si /-tpi / --Ri/2 fc 
JTQ I I 


dv 

~dp 

dv 

dp 


(s,<p) 


(s,<p) 


s ds dtp 


s ds dip 


vi rp(oo) ds 


rv 1 r 
JVo JPk 


dtp 


|Vv(j/)| d y, 


<Pi - V^o JC(k) 

where C(k) = {(p, <p), 2 ^ < P < ^t, <Pi < |<p| < <Po}- Hence, 

l06a - > 

IC(k) 


J 2 ^ 22fc / / |Vu(p)| 2 dpdx, 

<Pi - 7b “ iyfc ic(fe) 


and there is a constant C 3 > 0 independent of v such that 

■h < c 3 W |Vn(y )| 2 dp. 

fe^i^ c ( fc ) 

Note that every point (p, <p) in C(fc) lies in at most log 2 7 ^ sets C(2), l > 1. Therefore, 

J 2 < c 3 log 2 —!- [ \\7v{y)\ 2 dy ^ c 3 log 2 —!-V [ |Vu | 2 dx. 

-O-0 Jc 0-0 “ Jyfc 


(58) 


By the inequalities (l57l) and (l58l) . I 2 < C 4 X)fc>j /yfc |Vn | 2 dx for some constant C 4 > 0 independent of u. 
Hence, 

72 < 7EE» P_ ‘ +1 e l (r Y v>2dx ~ 

k^lp^k T£Y. k nA p }t(Y ’ 


(59) 
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Indeed, observe that the terms in the sum of (l59l) for which p = k are exactly k / rfc |Vu| dx. 

Therefore, (1521) and (l5!?l) yield, which C!)l) . which achieves the proof. □ 

Proof of Theorem 0 In the proof we will write < when there may arise in the inequality a constant that 
does not depend on the function v £ if 1 (flint) we consider. 

By Lemma[3]and Theorem fTOl for every v £ such that U|r £ £f 1 / 2 ( T), there exists Av £ H 1 { 

such that (Ari)|c> = v and 

v a 


IIAt 


inqR 2 ) 


< 


\m(c) 


LHC) 




Define C = Int(C' Uh ) (see Figure [I]). Introducing a cut-off function with support in the main hole T 
and using the operator A and Remark 1161 we can construct a linear extension operator Tq from iJ 1 (C') 
to 7J 1 (C' U T) such that T\ = 1 and for all v £ IF 1 (C) 

2 


L 


IVCFotOr dx 


< 


|Vw | 2 dx 


which also implies 


ll^oxll 


if 1 (CUT) 


< 


<c 


IHI 


dx, 


HHC) 


(60) 


(61) 


LHC) 


We will now define an extension v £ IF 1 (flint) of a function v as in Theorem [9l where flint is the convex 
hull of the domain flint- Recall that T is the main hole and {fa(T), a £ A} is the collection of the holes 
of the domain flint (see Figure 2|). Introduce the function v defined in flint by 

v:=v in flint, 

v \= T 0 (v o f a ) o / CT -1 in faiT), a £ A. 

By Lemma [3] and (1601) , we get the estimate 


/ |V(Jox)| 2 dx < J^ 2 k f 
JT k>1 J 7 


'-(^/dM + EE* 

k^Op^k 


p—k+1 




|Vu | 2 dx. (62) 


Indeed, k f v0 |Vu | 2 dx ( resp. k f c |Vu | 2 dx) is bounded from above by the terms for which k = p = 0 
(resp . k ^ 1 , p = k) in the second sum in (1621) . 

Observe that for every integer k ^ 1, C T' 7 where a k = (1,2,..., 2) e Ak- 2 , (recall that the sets T a 
have been introduced in Remark [3]). Therefore, 


f \v~{v) \ 2 dp < [ \v - {v) flrk \ 2 dx, 

Jik J 


(63) 


for i = 1 , 2 , where the constant in the inequality does not depend on k. Take a £ A n , one has 


/ 


U(T) 


|Vu | 2 dx = / |V(J- 0 (xo / ff ))| 2 dx 


< 


E* 


k-\-n 


\v - (v) 


SAL” 


dp 


k> 1 


SAL” k 

.p—k+1 f 

k^Op^k T£Z k nA p ■'ftT-P ’°) 


+ EE k 


|Vu | 2 dx, 


where we applied (RT^l) to the function v o f a , and we used (IMl) . The constant in the inequality does not 
depend on n. The notation or for r £ Ak stands for (er(l),..., cr(n), r(l),..., r(fc)) £ A n +k- 
We can write 


I |VT?| 2 dx = [ |Vu| 2 dx+^ [ |V5| 

J Lli n t -'Hint crGA StriT) 


dx 


< 


|Vu| dx + S 1 +S 2 , 
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where 


■si = EEE^/ lf J-(V) i 2 - 1 "- 

n^OaGAnk^l J ) 

s 2 = EEE « p_fc+1 E / o M 2 dx. 

crGAk^Op^k TGS k nA p 

We first deal with Si. Take er, r £ *4 and k,l ^ 1. Note that if erer fc = red, then fc = l and a = r. 
Therefore, 

Si < E 2 " E L (^)r-l 2 - H w |rllffi/a(r) ( 64 ) 

by (El). 

We are left with dealing with S 2 . Assume that 77 £ An and 77 = err with r £ E k (1 A p , then p ^ TV. Since 
the sets £ fe , fc ^ 1, are pairwise disjoint, this means that the term f^ |Vi>| da; appears at most N 

times in the sum S 2 . 

Moreover, we observe that p — k + 1 £ [1, N + 1], It can be seen that there is at most one quadruplet 
(cr', r', l, q) with o' £ A, t’ £ T, 1 f~l A q , l 1 and q ^ l distinct from (cr, r, k,p) such that 77 = o't' and 
p — k + \ = q — l + 1. As a consequence, 


JV +1 


< 2 £ £ £*”7 

N^0r]GA N m=l J U( Y °) 

2 




1 — ft 


E f |Vw| 2 da; = 


|Vw| 2 da; 
2 


1 — ft 


|Vu| 2 da;. 


^in 


Therefore, (15T1) and (1551) give 


INI/ri(n lnt ) ~ IMIff^nmt) + ll’ u |rll ff i/2 (r) - 


(65) 


Since flint is a polygonal domain, we can further extend v into a function Tv in 7J 1 (Z)), where J 7 is a 
linear operator satisfying C3l) . □ 


Proof of Theorem [8] We will prove separately points [I] and [777] in Definition [2] 

Proof of point\T\ Suppose that (u n ) weakly converges to u in L 2 (D). Without loss of generality, one can 
suppose lim a n (u n , u n ) is finite. We may further assume that there exists a subsequence, still called (u„), 
such that a n (u n , u n ) converges to some real number as 7i —> 00 ; as a consequence, there exists a constant 
c independent of n such that 

a n (u n , u n ) ^ c. (66) 

In particular, for all n, u n £ V n , which implies that u n £ V. Then, (l66l) implies that (un|n int ) is bounded 
in and that yJv n \Ju n iQ ext is bounded in L 2 (fl ext ). Therefore, there exists a subsequence that 

we still denote (u n ) such that 

• converges to U|n int weakly in H l (Pl int ), and strongly in L 2 (fl int ) 

• \fvn Vu ra |n ext converges weakly in L 2 (f2 ext ) to Vrt|Q ext (recall that v n converges to 1 almost every¬ 
where, so the weak limit of ypTf Vu ra |n ext must be Vzt|n ext .) 
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Moreover, exactly as in the proof of Theorem [6j we see that ^ry / r „ u 2 da: — > f r u 2 dp, as n —> 00 . We 
have proved point [T] 

Proof of point Take u £ L 2 (D). By (l24l) . we may assume that u £ V. We must construct (u n ) 
converging strongly in L 2 (D) such that (1271) holds. 

Recall that H is the height of the triangle T and that T° is a segment of the line { x 2 = 0}, T is a 
segment of the line { x 2 = 2 H}. We then introduce a sequence of smooth cut-off functions Xn(x 2 ) such 
that Xn{x 2 ) = 0 if x 2 <2 H - ^Irr, Xn{x 2 ) = 1 if x 2 £ [2 H - §c,2H] and that + ||Xn||L~ is 

bounded by a constant independent of n. 

We now define the functions u n by 


u n (x) = < 


Uint(x) 

ttext (x) 

£ ) (x) 


Xn(x 2 )£(Uj nt )(x) + (1 - Xn(x 2 ))Uext(x) 


Wx £ flint, 

Vx £ U, 

U U /-( T )> 

m>n cr^A-m 
n— 1 

Vx£ U U -M T )> 

m —0 cr^Am 


(67) 


where £ is the extension operator introduced in Theorem [9l u; n t = it|n int and u ex t = M|n ex f It is easy to 
check that u n belongs to the space V n and that the sequence ( u n ) strongly converges to u in L 2 (D). We 
claim that 


E E 

m —0 cr^Am 



V n {x) IV {Xn(£(uint) 


U ex t))| 2 dx ->■ 0, 


EE/ |V (^(lijnt) n e xt) | dx ^ 0. 

m>n creAm 7ct(X) 


( 68 ) 

(69) 


Indeed, we readily obtain (1S51) from the fact that the measure of U m > n f<r(T) tends to zero 

and the fact that ||f(uint)||ffi(D) and Hwextllt/qD) are finite. We obtain (1HS1) because ||v^YVxnlU^ is 
bounded uniformly with respect to n and Xn is supported in a region with vanishing measure, and because 
\\£{,u int )\\m(D) and ||uext||ffi(n ext ) are finite. 

Therefore, 


lim 

n—y oo 

I 




(|V< t | 2 + /3|<nt| 2 ) dx + [ K|V< xt | 2 + /3|< xt | 2 ) dx) 

t ^ ^ext / 


(|VM int | 2 + Puf nt ) dx + / (|V«ext| 2 + ^ext) d X, 


where u™ nt := w n |nr* t and M™ xt := u n |nj xt - Finally, from Proposition [T] jj4| f pn u 2 dx —> / r u 2 dp as 
n oo. Collecting all the above results, we obtain that lima rl (M r i,u rl ) = a(u,u) as n —> oo, thus point 
ii. 
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